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Quantum computation and communication rely 
on the ability to manipulate quantum states ro- 
bustly and with high fidelity. Thus, some form of 
error correction is needed to protect fragile quan- 
tum superposition states from corruption by so- 
called decoherence noise. Indeed, the discovery 
of quantum error correction (QEC)ii^ turned the 
field of quantum information from an academic 
curiosity into a developing technology. Here we 
present a continuous-variable experimental im- 
plementation of a QEC code, based upon entan- 
glement among 9 optical beams^. In principle, 
this 9-wavepacket adaptation of Shor's original 9- 
qubit scheme^ allows for full quantum error cor- 
rection against an arbitrary single-beam (single- 
party) error. 

QEC protocols eliminate uncontrolled errors that af- 
fect fragile quantum superposition states by encoding 
these quantum states into a larger, multi-partite entan- 
gled system. Errors occurring on a limited number of 
parties will leave the entanglement intact and so the orig- 
inal state may be retrieved by error syndrome recogni- 
tion followed by recovery operations. Shor proposed a 
concatenated quantum code to protect against arbitrary 
single-qubit errors, by encoding an arbitrary single-qubit 
state — a|0) + /3|1) into nine physical qubits 



encode 
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with |±) = (|0,0,0)± |1,1,1))/V2. Though reminis- 
cent of the redundant encoding in classical error correc- 
tion, the quantum code exhibits some clearly nonclas- 
sical features of which the most significant is the pres- 
ence of multi-party entanglement. The concatenation of 
three-party entangled states (|±)) into nine-party states 
enables one to correct both bit-flip and phase-flip er- 
rors. The latter type of error occurs only in nonclassical 
states. Remarkably, suitable error syndrome measure- 
ments would collapse an arbitrary error (including coher- 
ent superpositions of bit-flip and phase-flip errors) into 
the discrete set of only bit-flip and/or phase- flip errors. 
These discrete (Pauli) errors can be easily reversed to 
recover the original state. 



The continuous-variable version of Shor's 9-qubit 
codeii^ is the only code to date which can be determin- 
istically (unconditionally) implemented using only linear 
optics and sources of entanglement. Indeed, previous im- 
plementations of QEC were based on qubit codes, either 
in liquid-state NMR or linear ion trap hardware con- 
figurations. The liquid-state NMR experiments imple- 
mented QEC codes with up to five physical qubits^'^i^'i 
and in the ion trap experiment, a three-qubit code was 
realized^. Both configurations rely on nonlinear qubit- 
qubit coupling (in the form of nearest-neighbor cou- 
plings for NMR or via the collective vibrational mode 
for ion traps). Our experiment is the first implemen- 
tation of a Shor-type code, as the preparation of nine- 
party entanglement is still beyond the scope of exist- 
ing non-optical approaches and single-photon-based, op- 
tical schemes. Here, continuous- variable QEC— is real- 
ized using squeezed states of light and networks of beam 
splitters'^. Even this optical approach requires an opti- 
cal network three times the size as that used in earlier 
experimentsii to achieve the large-scale multi-partite en- 
tanglement for a 9-wavepacket code. 

In our scheme, as for the simplest QEC codes (whether 
for qubits or for continuous variables), a single, arbi- 
trary error can be corrected. Such schemes typically as- 
sume errors occur stochastically and therefore rely on 
the low frequency of multiple errors. Stochastic error 
models may describe, e.g., stochastic, depolarizing chan- 
nels for qubits, or in the continuous-variable regim e^'^'^'^ , 
free-space channels with atmospheric fiuctuations caus- 
ing beam jitter, as considered recently for various non- 
deterministic distillation protocols^^i^iii^ii^. For the 
continuous-variable QEC protocols, as realized in the 
present work, this type of error may be be suppressed 
in a deterministic fashion (see appendix F). The over- 
all performance of this family of QEC codes is then only 
limited by the accuracy with which ancilla state prepara- 
tion, encoding and decoding circuits, and syndrome ex- 
traction and recovery operations can be achieved. In the 
continuous-variable scheme, all these ingredients can be 
highly efficiently implemented; the finite squeezing of the 
auxiliary modes being the only limitation. This ancilla 
squeezing is linked with the presence of entanglement and 
it also determines whether the transfer fidelities exceed 
those of classical error correction (see appendices E and 

F). _ 

We begin with a description of the scheme in the limit 
of infinite squeezing, where the position x and momen- 
tum p of a harmonic oscillator (corresponding to a single 
optical mode of the light field) serve as the conjugate pair 
of observables used for the encoding 
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Through this 9-wavepacket code an arbitrary single-mode 
state IV-") = / dx^{x)\x) is encoded into nine optical 
modes. This perfectly encoded state is obtained by using 
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FIG. 1: 9-wavepacket quantum error correction cod&- for cor- 
recting an arbitrary error occurring in any one of the nine 
channels. The gray dotted hues represent the classical in- 
formation that is used to compute the necessary syndrome 
recovery operations. 



eight infinitely squeezed ancilla states. Finite squeezing 
of the ancillae leads to an approximate encoding, and 
hence lowers the fidelity of the QEC. 

Fig. 1 shows a schematic of our realization of the 9- 
wavepacket code. In the encoding stage, an input state is 
entangled with eight squeezed ancillae, each correspond- 
ing to an approximate '0' ("blank") state. After an error 
is introduced, the states are decoded simply by inverting 
the encoding. The eight ancilla modes are then measured 
(with a;-quadrature measurements performed in detectors 
1 and 4 and p-quadrature measurements in six other de- 
tectors), and the results of the measurement are used 
for error syndrome recognition. More precisely, these are 
the results of homodyne detection applied to the ancilla 
modes along their initial squeezing direction. 

The encoding stage consists of two steps in order to 
realize the concatenation of position and momentum 
codes^. First, position-encoding is achieved via a tritter 
T'in,ani,an47 that is two beam splitters (blue and green in 
Fig. 1) acting upon the input mode and two x-squeezed 
ancilla modes (anl and an4 in Fig. 1). The second step 
provides the momentum-encoding via three more fritters, 
with six additional p-squeezed ancilla modes (an2, an3, 
an5, an6, an? and an8 in Fig. 1). The overall encoding 
circuit becomes^S 



an7,aii8-^anl,an5,an6-' in,an2,an3-' in,anl,aii4 



(3) 



As the decoding stage merely inverts the encoding, the 
eight ancilla modes will remain all '0' in the absence of 
errors. In the presence of an error in any one of the nine 
channels, the measurement results of the decoded ancil- 
lae will lead to non-zero components, containing sufficient 
information for identifying and hence correcting the er- 
ror (see appendices for derivations and Table U for an 
error-syndrome map) . Similar to the qubit QEC scheme, 
where the conditional state after the syndrome measure- 



ments becomes the original input state up to some dis- 
crete Pauli errors, our conditional state coincides with 
the input state up to some simple phase-space displace- 
ments. Thus, it remains only to apply the appropriate 
(inverse) displacement operations in order to correct the 
errors. 

The detailed experimental setup for our 9-wavepacket 
QEC scheme is shown in Fig. 2. Eight squeezed vacua 
are created by four optical parametric oscillators (OPOs), 
which have two counter-propagating modes; thus, every 
OPO creates two individual squeezed vacua. The squeez- 
ing level of each single-mode squeezed vacuum state cor- 
responds to roughly 1 dB below shot noise. For pumping 
the OPOs, the 2nd harmonic of a cw Ti:Sapphire laser 
output is used. The syndrome measurements are per- 
formed via homodyne detection with near-unit efficiency. 

To apply a single error, a coherent modulation is first 
generated in a so-called error beam using an electro-optic 
modulator (EOM) ("modulated mode"). This beam 
is then superimposed onto the selected mode or chan- 
nel ("target mode") through a high-reflectivity beam 
splitter-1 with independently swept phase, resulting in 
a quasi-random displacement error. The error-correcting 
displacement operations (as determined by decoding and 
measurement) are then performed similarly, via an EOM 
and a high-reflectivity beam splitter, but now with phase 
locking between the modulated and target modes along 
either the x or p axis, as appropriate. 



TABLE I: Error syndrome measurements. LO phase: quadra- 
ture at which the local oscillator phase of the homodyne de- 
tector is locked, ES: equal signs, DS: different signs. 
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FIG. 2: Experimental setup of the 9-wavepacket quantum error correction; PBS: polarization beam splitter, PPKTP: periodi- 
cally poled KTiOP04, HBS: half (symmetric) beam splitter, HWP: half wave plate, ND: neutral density filter, PZT: piezoelec- 
tric transducer, BHD: balanced homodyning, SHD: self-homodyning, OPO: optical parametric oscillator, MCC: mode-cleaning 
cavity, LO: local oscillator, ISO: optical isolator, EOM: electro-optic modulator. 



Fig. 3 shows some examples for error syndrome mea- 
surement results. Here, the input state is chosen to be 
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FIG. 3: Error syndrome measurement results. (I) A random 
displacement error is imposed on channel 1. (II) A random 
displacement error is imposed on channel 9. A two-channel 
oscilloscope is used measuring the outputs of detectors 1 and 

4, 2 and 3, 5 and 6 and 7 and 8. (a) output signal of detector 
1, (b) detector 2, (c) detector 3, (d) detector 4, (e) detector 

5, (f) detector 6, (g) detector 7, (h) detector 8. 



a vacuum state. A random displacement error in phase 
space is imposed on channel 1 (Fig. 3(1)) and on channel 
9 (Fig. 3(11)). A two-channel oscilloscope is used to mea- 
sure the outputs of pairs of detectors (1, 4), (2, 3), (5, 
6), and (7, 8). Comparing the results of Fig. 3(1) with 
Table 1, one can identify an error occurring in channel 
1, since only detectors 1 and 2 have non-zero outputs. 
The outputs from detectors 1 and 2 correspond to the 
desired x and p displacements, respectively. Similarly, 
from Fig. 3(11), we can recognize that an error has oc- 
curred in channel 9. Here, detectors I, 4, 7, and 8 have 
non-zero outputs and the outputs of detectors 1 and 4, 
as well as 7 and 8 have equal signs (distinguishing it from 
the case of an error in channel 8, for which outcomes 7 
and 8 have different signs). 

Fig. 4 shows two examples of QEC results, comparing 
output states with and without error correction, and with 
and without squeezing of the ancilla modes. In Fig. 4(1), 
an error was introduced in channel 1. The local oscillator 
(LO) phase of the homodyne detector was tuned to detect 
the X quadrature of channel 1. Similarly, in Fig. 4(11), 
the error was introduced in channel 9 and the LO phase 
is locked to the p quadrature. For ease of experimental 
implementation, only the measurement outcomes of de- 
tectors 4 and 8 were fed forward to the error correction 
step. In principle, using the combined outputs of detec- 
tors 1 and 4 for x and detectors 7 and 8 for p would yield 
even higher fidelities. 

The quality of the error correction can be assessed via 
the fidelity F = (i/'inl Pout | "0111), where j-^in) represents the 
input state and pout corresponds to the output state of 
the error correction circuit^liSSiiS^. Here the fidelity is 
calculated as 



F = 



V(l + 4((Aiout)2))(l + 4((Apout)2)) 
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FIG. 4: Results of quantum error correction. (I) A random 
phase-space displacement error is imposed on channel 1. The 
LO phase of the homodyne detector is locked to the x quadra- 
ture. (II) A random displacement error is imposed on channel 
9. The LO phase of the homodyne detector is locked to the 
p quadrature. In each case, four traces are shown compari- 
son: (i) Homodyne detector output without error correction 
(no feed forward step), (ii) Error correction output without 
squeezing, (iii) Error correction output with squeezing, (iv) 
shot noise level, (a) Single scan of a spectrum analyzer with 
zero span mode. 2 MHz center frequency, 30 kHz resolution 
band width and 300 Hz video band width, (b) 30 times aver- 
age of traces (ii-iv) above. 



where Xout and Pout are quadrature operators of the out- 
put field. For example, in the case of an error in channel 
1, the output quadrature operators become 



1 



1 



_ _ -(0) -r2 

Pout — Pin ^Pan2^ 



(5) 



where iin, Pin, ^ini: and p)^^2 are quadrature operators 
of the input field and the ancilla vacuum modes, and 
ri are squeezing parameters for ancilla i. In the ideal 
case of ri —^ oo, unit fidelity is obtained, with output 
states approaching the input states. For zero squeezing, 
Eq. (O yields an excess noise of i and i for the x and 
p quadratures, corresponding to 1.76 dB and 0.67 dB of 
output powers, respectively (see Table HI)) . 

Eq. (4) can be used to translate the measured noise 
level values from Table HTl into fidelity values. Indeed, for 
every possible error introduced (in any of the channels) 
the fidelity after error correction exceeds the maximum 
values achievable for the scheme in the absence of ancilla 
squeezing. For example, for an error in mode 1, a fidelity 
of 0.88±0.01 was achieved (exceeding the classical cutoff 
of 0.86). Similarly, for an error in channel 9, we obtain 



s(o) 



TABLE II: Output noise power of QEC circuit in dB, rela- 
tive to the shot noise level. Perfect error correction therefore 
corresponds to dB. SQV: squeezed vacua. 
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a fidelity of 0.86±0.01, exceeding a cutoff of 0.82. (The 
lower cutoff takes into consideration that only two of the 
four non-zero components are used.) The better-than- 
classical fidelities for errors in any one of the nine chan- 
nels are indirect evidence of entanglement-enhanced error 
correction (see appendices). By comparison, in complete 
absence of any error correction, i.e., without reversing 
displacement errors (including the zero-squeezing case; 
for an application of such "classical" error correction, see 
appendix F), fidelity values under 0.007 ±0.001 were ob- 
tained. 

In conclusion, we experimentally demonstrated a Shor- 
type quantum error correction scheme based upon entan- 
glement among nine optical beams. The entanglement 
is used for deterministically generating a concatenated 
code, allowing for the correction of arbitrary errors in 
any one of nine communication channels. In the experi- 
ment, evidence is obtained for an entanglement-enhanced 
correction of displacement errors; a further increase of 
the small enhancement of the current implementation 
would only require higher squeezing levels of the resource 
states. Our experiment represents the first demonstra- 
tion of quantum error correction beyond qubits (and 
specifically for continuous variables). The scheme may 
be useful for any application in which stochastic errors 
occur such as free-space communication with fiuctuating 
losses and beam pointing error o^^i^^'^^i^^ . The ability to 
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implement QEC in an optical network of this size rep- 
resents a significant step towards the manipulation and 
application of large-scale multi-partite entanglement for 
quantum information processing. 
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APPENDIX A: ENCODING 



Equation (3) describes a nine-port device acting upon 
the signal input mode, two x-squeezed ancilla modes 
("anl" and "an4" in Fig. 1 of main body), and six p- 
squeezed ancilla modes ("an2", "an3", "an5", "an6", 
"an7", and "an8" in Fig. 1). Labeling the nine input 
modes by subscripts one through nine, we obtain the out- 



put quadrature operators of the encoded state, 
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Note that with respect to these subscripts, eq. (3) can be 
expressed by T'789J456Ti23Ti47 for modes 1 (signal input), 
2 ("an2"), 3 ("an3"), 4 ("anl"), 5 ("an5"), 6 ("an6"), 7 
("an4"), 8 ("an7"), and 9 ("an8"). 

The encoded state exhibits the following quadrature 
quantum correlations in the case of nonzero squeezing, 
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In the limit ris 00, the quadrature operators become 
perfectly correlated. 
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These correlations are analogous expressions to the eight 
stabilizer conditions of the Shor qubit code (where for 
continuous variables, Pauli operators are replaced by 
Weyl-Heisenberg phase-space operators). Note that 
these correlations hold for any signal input state, i.e., 
for any resulting "code words" , again similar to the sta- 
bilizer conditions for qubits. In order to obtain a suffi- 
cient set of entanglement witnesses for verifying a fully 
inseparable nine-party state, additional quadrature cor- 
relations must be considered; these extra correlations are 
expressed in terms of the "logical" quadratures in the 
code space which depend also on the signal state (see 
appendix E). 



APPENDIX B: DECODING AND CORRECTION 

Random phase fluctuations are transferred onto one 
selected beam of the encoded state, leading to random 
phase-space displacements of one of the nine optical 
modes. This effect can be described by adding error 
quadrature operators to the corresponding mode fc, AfcX^ 
and Afep^, where the parameter will be set to one 
for the single mode of the noisy quantum channel and 
otherwise chosen to be zero. After the decoding step, 
^147^123^456-^7^91 Outgoing quadrature operators be- 
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Modes two through nine are measured via suitable homo- 
dyne detectors, i.e., the local oscillator phase is adjusted 
to detect those quadratures which are quiet if there was 
no error. After the corresponding feedforward operations 
on the first mode, the signal input state will be recov- 
ered in mode 1 up to the finite squeezing from the ancilla 
modes. 

For example, in the case of an error transferred onto 
mode 1, Afe = 4i, 

— .Till I 2*^1' 
1 

Pi = Pin + gPl, 

(B2) 

only for detectors 1 and 2 (see Fig. 2 of main body), 
measuring x'^ (position of "anl" ) and p'2 (momentum of 
"an2"), respectively, results clearly different from zero 
(coming from the error) are obtained. All the remaining 
detectors show results around zero. In order to correct 
the error, mode 1 is displaced according to 

1 

V2 
1 



Pi ^ Pi - ^P2: 



(B3) 



using eqs. (|B1[) . 

Similarly, in the case of an error transferred onto mode 
9, Afe = (5feg, we have 

■^1 ~ -^in + 3'^9' 
Pi = Pin + gPg- 

(B5) 

Now the only nonzero outputs occur at detectors 1 and 
4 (Fig. 2), measuring x'^ (position of "anl") and (po- 
sition of "an4"), respectively, and at detectors 7 and 8, 
measuring pg (momentum of "an7" ) and pg (momentum 
of "an8" ) , respectively. Possible correction displacements 
are 



x\ + 



-•7' 



(B6) 



for X, and 



P'l 



P'l 



Ip',, 

, V2 . 



Pi ^ Pi + — Pg> (B7) 
for p. These corrections result in the output quadratures 



•^out 
■^out 

Pout 

Pout 



x/2t^°^ 

V ^-^anl 



(0) ^-r. 



;(0) 



<^ ^ an4 ^ : 
Pin+\/-Pa„7e > 



' I ^ -(0) -rs 
Pin + g Pan8^ ' 



(B8) 



always nearly recovering the signal input state. 

Similar calculations yield the quadrature operators for 
the output state of mode 1 after the error correction pro- 
tocol in the case of an error on modes two through eight; 
for an error on mode 2, 



1 



'^out.detl 



Pout,dct2 — Pin 



^^out.detS 



^■^anl 

^2 „(o) 



%/2 



Pan3^ ' 



(B9) 
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for an error on mode 3, 

3^out,detl — ^in 

Pout,dot2 = Pin 

Pout, dots — Pin 

for an error on mode 4, 

•^outjdctl — -^in 

•^out,det4 — •^in 

Pout, dots = Pin 

for an error on mode 5, 

•^out,detl — -^in 

^out,dct4 — -^^in 

Pout, dots — Pin 

Pout,dot6 — Pin 

for an error on mode 6, 

^out,dotl — ^in 

^out.dot4 — ^in 



V2' 



.(0) 
^anl^ 



Pan2^ ! 



%/2 .(0) 

2 Pan3^ ' 



(0) g-.. 



■^an4"^ ' 
PanS'^ ' 



V6 



2 -(0) -rs 
gPanS^ 



V2 



JO) g-re 
Pan6^ ' 



Pout, dots — Pin + 

Pout, dote = Pin 

for an error on mode 7, 

■^out,dotl — -^in 

•^out,dot4 — -^in 

Pout, dot? = Pin 

for an error on mode 8, 

•^out,dotl — -^in 

*^out,det4 — -^in 

Pout, dot? = Pin - 

Pout, dots = Pin ' 



;(0) 



I JO) „~r, 
g-fanS^ ' 



V2 .(0) -rs 
2 Pan6^ ' 



^,(0) p-r-4 
■^an4"^ ' 



JO) p-r-r 
PanT^ ' 



x/6 



V24°n\e--, 



;(0) 



gPan?^ ' 



V2 



.(0) - 
Pan8^ 



(BIO) 



(Bll) 



(B12) 



(B13) 



(B14) 



The additional subscripts "detl", etc., indicate which 
detector outcomes are used for the correction displace- 
ments. These detectors (see Fig. 2) measure the quadra- 
tures x'^ ("detl"), p'^ ("det2"), p'^ ("detS"), x'^ ("det4"), 
p^ ("det5"),p^ ("det6"),p^ ("det7"), and ("detS"). 

Because of the freedom in choosing the correction dis- 
placements, there is always an optimal feedforward op- 
eration. For example, in the case of an error on mode 
2, 



•^out 
Pout,det2 
Pout,dct3 



1 

V2 



x(0) 



Pil 



JO) 

Pan2e , 



V2 .(0) 

Pin 2 PanS'S ' 



(B16) 



we obtain the following excess noise for the output state, 



((p, 
((p. 



'out,dct2 ) 



'out,dct3 ) 



However, for r2 = ^3 
ation leads to 



{{iin?) 
{{Pin?) 
{{Pin?) 



1 1 
— • — f 

2 4 
2 1 

2 1 



-2ri 



-2r2 



(B17) 
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r, the optimal feedforward oper- 



Pout,opt — Pin 

corresponding to 



1 -(0) -r _ V2 .(0) _ 



((p. 



out, opt ) 



{{Pin?) + \-\e 



2r 



(B18) 



(B19) 



(B15) 



which is the same as for the case of an error on mode 
1. For unequal squeezing, r2 ^ r^, the optimal feedfor- 
ward depends on the squeezing values. Therefore, in the 
current experiment, we use only the output of detector 3 
for the feedforward. Table IB] shows which outputs of the 
homodyne detectors are used for error correction. 



APPENDIX C: RESULTS OF ERROR 
SYNDROME MEASUREMENTS 

Fig. [5] shows error syndrome measurement results. 
Here, the input state is a vacuum state. This case is 
also described in the main body of the paper. A ran- 
dom displacement error in phase space is transferred onto 
quantum channels one through nine for A-I, respectively. 
With Table 1 in the main body of the paper and Fig. [5] 
here, we can decide which quantum channel is subject to 
an error and derive a corresponding feedforward opera- 
tion to correct the error. 



9 







1W f — 1 — 

u ^ 

■K 

-llB E — 1 — 




1W f — 1 — 
SO I'l 




■V) 




1«1 f — 1 — 

D 








i» f — 1 — 









IW F — 1 — 
M 




■V) 




H» f 1 

M 13) 

[1 




■V) 




1W1 f — 1 1 1 1 1 1 1 1 1 — □ 


• S(l 



:?\l'U\r\ 



I 

I 



I »['•' 



1 



i i 
i ™ 



JMVJW\ 



i to 



r 

iWr-" — 




» ^"^ 

i i 

■» ^ 




r 






VAAaA 










i»F— ■ 


1 




»|« 





-liM E — 1 — 1 1 — 1 — 1 — 1 








E — 1 — 1 — 1 — 1 — 1 — 1 








.100 E — 1 — 1 — 1 — 1 — 1 — 1 


,,, 


-HW t — 1 — 1 1 — 1 — 1 — 1 




.100 E — 1 — 1 — 1 — 1 — 1 — 1 








GU ' 




.100 E — 1 1 1 1 1 1 











I 4 



1 " 

^ IM 




1 " 




1. SI 

t: 


(If) ■ ' ' ' 
















s 




* z 








„ ,1. 

1. M 

f ^ 




I : 

S g 

i : 




? 

f ^ 

l« 


w ■ 


f ™ 










1. M 


(!> ' ' ' ' 




(Tl ' ■ ' ' 


^ -M 
■' iM 




1,: 




1 M 

f 


|!f) ■ ' ' ' 


i : 

i» 


[h5 ' ■ ' ' 


E Z 
i 




!^ 
1 -™ 





a n X M 



14 » H 
TbTMlHI 



FIG. 5: Results of error syndrome measurements. A-I correspond to the cases of an error in quantum channels 1-9, respectively. 
(a)-(h) correspond to outputs from homodyne detectors 1-8, respectively. 



APPENDIX D: RESULTS OF ERROR 
CORRECTION 

Fig.[S]shows the results of error correction. The results 
are summarized in Table 2 in the main body of the paper. 



APPENDIX E: THE ROLE OF MULTIPARTITE 
ENTANGLEMENT 

The encoded nine- mode state, as created in the current 
experiment and described by eqs. (jAip . approaches the 



following state in the limit of infinite squeezing, 

IV'encodc) = ^ I dp dpidp2dp3^{p)e~^'P^P'+P'+P''> 
X\P1,P1,P1,P2,P2,P2,P3,P3,P3)- (El) 

Clearly, even for infinite squeezing and perfect encod- 
ing, the inseparability properties of the total nine-party 
state depend on the signal input wave function 'il'{p). 
In particular, for ^^^{p) = 5{p), we obtain IV'encodc) — 
J dp\p,p,p)® J dp |p, p, p) (X) J dp \p, p, p) , which is clearly 
not fully nine-party entangled, but rather a product state 
of three fully tripartite entangled GHZ-type three-mode 
states. So in order to obtain full nine-party entangle- 
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mm 



FIG. 6: Results of error correction. A-I correspond to the cases of an error in quantum channels 1-9, respectively. The LO 
phase of the homodyne detector is locked at x or p, which is indicated after the capital letters. Trace numbers are the same as 
in Figs. (5) and (6) in the main body of the paper. 



ment, the input state should not correspond to an in- 
finitely x-squeezed state (corresponding, after Fourier 
transform, to 4'ip) = ^(p))- Similarly, an infinitely p- 
squeezed input state leads to vanishing GHZ-type corre- 
lations within each of the three triplets, but it has ex- 
cellent GHZ-type correlations between the three triplets. 
For an input state between these two extremes, for in- 
stance, a vacuum input state as used in the experiment, 
we obtain quadrature correlations of both types, poten- 
tially leading to full nine-party entanglement. 

In order to witness full nine-party entanglement, in ad- 
dition to the correlations of eqs. (|A3|) . p-correlations be- 
tween the triplets and x-correlations within each triplett 
are required. Equations (|A3p only describe ^-correlations 
between the triplets and p-correlations within each 



triplet. The missing correlations are of the type of 
Pi + Pa + P7 — > and xi + £2 + £3 0. These linear 
combinations correspond to the "logical" quadratures in 
the code space, 

X = xi + X2 + X3 = Xin + V^xfj^^e^""^ , 

P = P1+P4+P7 ^ Pin + ^pfn2P~''^ 

+ + /l^^i^e- , (E2) 

obviously depending on the signal input state. Only 
for infinite ancilla squeezing, the encoding is perfect, 
X — x-yn and P = pi^- The excess noise in each quadra- 
ture is 2 X e~^''/4 for equal squeezing of the ancilla 
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TABLE III: Homodyne detector outputs for feedforward in 
the current experiments. 



channel with 


quadrature 


detectors for 


an error 




feedforward 


1 


X 


1 




V 


2 


2 


X 


1 




P 


3 


3 


X 


1 




p 


3 


4 


X 


4 




p 


5 


5 


X 


4 




p 


6 


6 


X 


4 




p 


6 


7 


X 


4 




p 


7 


8 


X 


4 




p 


8 


9 


X 


4 




p 


8 



modes. In this case, an infinitely x-squeezed input state 
would lead to excellent intra-triplet a;-correlations; an in- 
finitely p-squeezed input state, favorable for good inter- 
triplet p-correlations, leads to vanishing intra-triplet x- 
correlations. With a vacuum input state, as used in the 
experiment, we have both types of quantum correlations 
for nonzero squeezing of the ancilla modes. Similar quan- 
tum correlations also exist for the combinations 



X4 -I- .T5 -I- Xe 
X7 + Xs + Xg 

P2+P5+ P8 



1 



2'^an4^ 



1 ^(0) . /3 .(0) -r4 



V2 

- _ 1 -(0) -r2 
Pin ^ygPan2^ 

1 

-JO) 



2'^an4^ ' 



JO) -r, 
2^^an3'^ 



2Pan6^ 
1.(0) 



P3+P6+P9 = 



1 -(0) -r. 
Pin - "^?'a„2e - \j T;Pln3e 

_J_^(0).-5_./i^(0) ^-re 



V6 V 2 

J_r5(o)p-7_,/i^(o) „ 
^Pan7e Y ^P^nae 



(E3) 

The total set of quadrature quantum correlations can 
be sufficient for a fully inseparable nine-party entangled 
state. The corresponding nine-party entanglement wit- 
nesses lead to the known criteria for multi-party insepa- 
rability of continuous- variable states^^. In order to verify 
three-party inseparability within each triplett P123, /5456, 
and P789, we have, for /5i23, 



([A(pi - 
([A(P2 - 

for P456, 

([A(P4 - 
([A(P5 

for (0789, 



P2)]2) + {[A{X1 +X2+ gia X3)]^) < 1, 
- P3)f) + ([A(.9ifc X1+X2+ X3)f) < 1, 



(E4) 



-P5)]'> 

-Pe)?) 



{[A{x4 + X5 
([A(526 £4 + ^5 



xe)?] 



< 1, 

< 1, 



(E5) 



mP7 ' 

([A(P8 



Ps)?) + {[A{X7 +XS+ gsa X9)?) < 1, 
- P9)]2) + ([A(53b £7 + ^8 + Xg)?) < 1. 



(E6) 



The "gains" gia, etc., can be used to optimize these 
conditions. In order to rule out a state of the form 

J-") ^ Ji) ^ Ji) = 



E^ VrP\i3 ® P\k ® P789 = E^ VrPa 

further criteria, for example, 



' pI'^ ® Pc^ , we need 



[A(pi + p4 + .911 pi)f) 
-l-([A(£i + gi2 X2 -I- gi3 x-i ■ 

[A(.g2iPi -^P^+Pi)f) 

-|-([A(£4 -f- 322 £5 + 323 ^6 ■ 



x^ - gi4 X5 - gi5 xq) 



X7 



.924 2^8 - ff25 2^9 ) 

which describe the inter-triplet correlations. 



(E7) 

< 1 

(E8) 

< 1 

Equa- 



tion IIE7I) rules out the forms J^i ViPa ^ ® Pbl ^^'^ 



(i) 



' Pac ; eq. 



(|E8|) rules out the forms ijipl:^ 

p'^l and J^iViPb^^ ® Pac- Thus, any form of separability 
between the triplets a, 6, and c can be ruled out. The 
inter-triplet conditions can be understood as GHZ-type 
correlations of modes 1, 4, and 7 after LOCC operations; 
namely, x-measurements of modes 2, 3, 5, 6, 8, 9 and the 
corresponding displacements of modes 1, 4, and 7. 

In the experiment, it was verified that in any of the 
nine cases of an error in any one of the nine channels, 
the classical cutoff (zero-squeezing limit) was exceeded. 
This confirms that all 8 ancilla modes are in a squeezed 
state (see Table 1 and Fig. 2), as the quadrature noise 
of every ancilla mode contributes to the excess noise of 
the corrected signal for some of the detector results used 
for feedforward. This squeezing translates into nonclas- 
sical correlations for all combinations in eqs. (|A2|) . (|E2p . 
and (|E3p (with a vacuum input state). The set of quadra- 
ture combinations corresponds to the "unit-gain" version 
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of the entanglement witnesses in eqs. (jE4l) . (jE5|l . (|E6p . 
(|E7[) . and (|E8p . In order to satisfy the witness inequah- 
ties, in particular, for small squeezing values (as those of 
roughly 1 dB in the experiment), non-unit gain must be 
chosen. Although these non-unit gain combinations have 
not been measured directly in the quantum error correc- 
tion experiment, the nonclassicality in all the unit-gain 
combinations may be interpreted as an indirect confir- 
mation of the presence of nine-party entanglement. 



The sign of the displacement error is fixed and known, 
e.g., without loss of generality, X2 > 0. Note that more 
general errors, including non-Gaussian x-errors, could be 
considered as well. 

Now in order to encode the input state, we use two 
ancilla modes, each in a single-mode a;-squeezed vacuum 
state, represented by 

W^^,{xk,Pk) = - exph2e+2''4 - 2e"2>2] ^ (^4) 



APPENDIX F: APPLICABILITY OF 
CONTINUOUS- VARIABLE CODES 

The continuous-variable nine-mode code corrects an 
arbitrary error occurring in any one of the nine channels. 
Similar to the qubit case, for realistic scenarios, we should 
consider imperfect transmissions in every channel under 
the reasonable assumption of errors acting independently 
in all the channels. The performance of the code can then 
be evaluated by comparing the transfer fidelities for the 
encoded scheme with a direct transmission of the signal 
state through a single noisy channel^^. 

Using the example of a 3-wavepacket code, we shall 
demonstrate that for certain stochastic error models, 
the continuous-variable code leads to a dramatic im- 
provement of fidelity even when the errors occur in 
every channel"'^'^. In this case, the errors should cor- 
respond to a;-displacements or any errors decompos- 
able into a;-displacements (including non-Gaussian "x- 
errors"). A code for correcting arbitrary errors includ- 
ing non-commuting x and p-errors is obtainable, for in- 
stance, by concatenating the 3-mode code into a 9-mode 
code, as implemented in the current experiment. The ap- 
propriate error models are reminiscent of the most typ- 
ical qubit channels such as bit-fiip and phase-fiip chan- 
nels. In the continuous-variable regime, these types of 
stochastic errors would map a Gaussian signal state into 
a non-Gaussian state represented by a discrete, incoher- 
ent mixture of the input state with a Gaussian (or even 
a non-Gaussian) state. 



(Fl) 



Here, the input state described by the Wigner function 
Win is transformed into a new state Wenor with prob- 
ability 7; it remains unchanged with probability 1 — 7. 
A special case of the above channel model is an erasure 
channel^. The generalized erasure model here may find 
applications in free-space communication with fiuctuat- 
ing losses and beam point jitter effect a-^^'^''i^^ . 

As an example, we will consider a coherent-state input, 
\ai) = \xi + ipi), described by the Wigner function, 

Winixupi) ^-exp[~2{xi~xif~2{pi-Pif].{F2) 

Moreover, we assume that the effect of the error is just 
an a:-displacement by X2 such that 



with squeezing parameter r and k = 2,3. The total three- 
mode state before encoding is 

W{ai,a2,a3) = T4^in(a;i,Pl)Wanc(2:2,P2)W^anc(a:^3,P3), 

(F5) 

with aj ~ Xj + ipj, j = 1,2,3. The encoding may 
be achieved by applying a "fritter", i.e., a sequence of 
two beam splitters with transmittances 1 : 2 and 1:1. 
The total, encoded state will be an entangled three-mode 
Gaussian state with Wigner function. 



VFcnc(ai, ^2, as) = 



(F6) 



X exp \ — 2 



—^[Xi+X2+ X3 



Xi 



-—e 
3 



1 / 
73 

,+2r 



(pi - P2) + {P2 - Ps) + [pi - PaY 
2 



[Pl +P2+P3 



Pi 



[xi - X2f' + {X2 - Xzf + {Xi 



X3, 



A) 



Worror(a;i,Pl) = Win(xi - X2,Pl) . 



(F3) 



Now we send the three modes through individual chan- 
nels where each channel acts independently upon every 
mode as described by Eq. (jFl[) with Wcrror corresponding 
to an a;-displacement by X2. As a result, the three noisy 
channels will turn the encoded state into the following 
three-mode state, 

W^cnc(ai>"2,a3) (F7) 

= (1 - 7)^Wcnc(ai,a2,a3) 

+7(1 - lYWc^c{xi - X2+ ipi, q;2, as) 

+7(1 - 7)^M/onc(ai, a;2 ~ X2 + ip2, "3) 

+7(1 - jYWenciai,a2,X3 - X2 + ips.) 

-1-7^(1 - 7)We„c(xi -X2+ ipi,X2 - X2 + ip2,a?,) 

-1-7^(1 - 7)Wo„c(a;i -X2+ ipi, 02^X3 - X2 + ipz) 

-1-7^(1 - 7)Wo„c(ai,a;2 - X2 + ip2,X3 - X2 + ipa) 

+l^Wcnc{xi -X2+ ipi,X2 - X2 + ip2, X3 - X2 + ips) ■ 

Note that we assumed the same x-displacements in every 
channel. 

The decoding procedure now simply means inverting 



13 



the tritter, which results in 

Wdoc(ai,a2,a3) (F8) 

1 



+7(1 - 7) Win ( Xi - —f=X2,Pl 



xWanc - ^ ■^X2,P2 j Wane (2:3,^3) 

+7(1 - -ffWin (^Xi - ^*2,Pl 

xWanc ^2^2 + "^^2,^2^ W^anc(a;3 " -^52,^3) 
+7(1 - jfWin (^Xi - ■^X2,Pl 

xWanc ^2^2 + -^^2,^2^ Wi,nc{x3 + -^X2,P3) 



+7 (1 - iWin \^Xi - —^X2,Pl 

xWanc ^2^2 " ^^^2,^2^ W^anc(2;3 " ■^X2,P3) 
+7^(1 - j)Win (^Xi - -^S2,Pl 

xWanc ^2^2 ^ -^X2,P2^ W^anc(2;3 + -^X2,P3) 



+7^(1 - 7)W^i„ ( Xi - -^X2,Pi 



xWanc ^X2 + ^ -X2,P2 j VFanc (2:3 , P3) 

+7^T^in (2:1 - VSX2,P1^ 

X Wane {X2,P2) Wa„c(x3,P3). 

By looking at this state, we can easily see that 2;- 
homodyne detections of the ancilla modes 2 and 3 
(the syndrome measurements) will almost unambigu- 
ously identify in which channel a displacement error oc- 
curred and how many modes were subject to a displace- 
ment error. The only ambiguity comes from the case 
of an error occurring in every channel at the same time 
(with probability 7'^), which is indistinguishable from the 
case where no error at all happens. In both cases, the two 
ancilla modes are transformed via decoding back into the 
two initial single-mode squeezed vacuum states. All the 
other cases, however, can be identified, provided the ini- 
tial squeezing r is sufficiently large such that the displace- 
ments cx X2, originating from the errors, can be resolved 
in the ancilla states. 

The recovery operation, i.e., the final phase-space dis- 
placement of mode 1 depends on the syndrome mea- 
surement results for modes 2 and 3 which are consistent 
with either undisplaced squeezed vacuum states ('0') or 
squeezed vacua displaced in either or '— ' x-direction. 



The syndrome results for modes 2 and 3 corresponding 
to the eight possibilities for the errors occurring in the 
three channels are (0,0) for no error at all, (+,0) for an 
error in channel 1, (— ,-|-) for an error in channel 2, (— ,— ) 
for an error in channel 3, (-l-,-l-) for errors in channels 1 
and 2, (+,— ) for errors in channels 1 and 3, (—,0) for 
errors in channels 2 and 3, and, again, (0,0) for errors 
occurring in all three channels. 

In the limit of infinite squeezing of the ancilla modes, 
the ensemble output state of mode 1 (upon averaging over 
all syndrome measurement results X2 and 2:3 including 
suitable feedforward operations) can be described as 

(1 - 7')Wi„(2;i,pi) + 7'Win (2:1 - V32;2,pi) . (F9) 

This output state emerges, because in almost all cases, 
the feedforward operations turn mode 1 back into the 
initial state (in the case of finite squeezing, only up to 
some Gaussian-distributed excess noise depending on the 
degree of squeezing used for the encoding). The only case 
for which no correction occurs is when errors appear in 
every channel at the same time, at a probability of 7'^. 
In this case, the initial state remains uncorrected, with 
an x-displacement error of y/^X2- 

We see that a fidelity of 1 — 7'^ can be achieved, assum- 
ing a;2 ^ 1 (for smaller X2, the fidelity would even exceed 
1 — 7^, but those smaller X2 may be too hard to detect 
at the syndrome extraction, depending on the degree of 
squeezing, see below). Note that this result implies that 
the encoded scheme performs better than the unencoded 
scheme (direct transmission with indirect = 1 ~ 7) for any 
< 7 < 1. In other words, by employing the quantum 
error correction protocol, the error probability can be 
reduced from 7 to 7"^. The continuous- variable scheme 
in this model is more efficient than the analogous qubit 
repetition code and it does not require error probabilities 
7 < 1/2 as for the case of qubit bit-flip errorsi^. 

We may now consider two different regimes for the 
error displacements X2- First, the regime e~^'"/4 < 
2;2 < 1/4, corresponding to small displacements below 
the shot noise limit; these can only be resolved provided 
the squeezing is large enough. In the limit of infinite 
squeezing r — s- 00, arbitrarily small shifts can be detected 
and perfectly corrected (with zero excess noise in the out- 
put states corresponding to unit fidelity). 

Secondly, the regime ^2 » 1. For these large shifts, 
even zero squeezing in the ancilla modes (i.e., vacuum 
ancilla states) is sufficient for error identification. Even 
with r = 0, the syndrome measurements still provide 
enough information on the location of the error and, to 
some extent, on the size of the error. We may refer to 
this kind of scheme as classical error correction (CEC), 
corresponding to the "classical cutoff" used as a classical 
boundary in the main body of the paper. This classical 
cutoff depends on the particular encoding and decoding 
circuit used; in the experiment, it is the same circuit as 
that employed for quantum error correction (neither of 
these are necessarily optimal). 
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CEC for large shifts .f 2 ^ 1 (the regime of the ex- 
periment) works fairly well. In fact, the fidelity values 
without CEC drop to near-zero fidelities, as measured 
in the experiment, F < 0.007 ± 0.001. Experimentally, 
this CEC is a highly nontrivial task and it is needed to 
achieve reasonable transfer fidelities. Nonetheless, the 
CEC scheme results in excess noise for the output state 
coming from the feedforward operations based on the 
fluctuating syndrome measurement results. By employ- 
ing squeezed-state ancilla modes, this excess noise can 



be reduced (down to zero for infinite squeezing). In this 
case, the scheme operates in the quantum regime. Sig- 
nificantly, in the experiment, non-commuting errors have 
been corrected, which means that CEC will always result 
in some excess noise. Although the margin of the demon- 
strated quantum error correction (on top of the CEC) is 
rather small, the experimental data provide clear evi- 
dence that CEC has been outperformed by the quantum 
scheme. 



